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ABSTRACT 

A well-known  result  of  S.  N.  Bernstein  states  that  a Chebyshev 


quadrature  formula  of  the  form 

/ 1 l «V  + Vf>’  \ real  ' 

-1  k=l 

can  have  algebraic  degree  of  exactness  p = n only  if  1 < n < 7 or 
n = 9 . The  nodes  t^  are  necessarily  symmetric  with  respect  to  the 
origin,  so  that  in  fact  p = 2[n/2]  + I . If  symmetry  of  the  nodes  is  im- 
posed, it  is  known  from  work  of  Gautschi  and  Yanagiwara  and  others  that 
next-to-highest  algebraic  degree  p = 2[n/2]  - 1,  beyond  the  classical 
cases  above,  can  be  attained  only  when  n = 8,  10,  11,  and  13.  For 
these  values  of  n,  '‘optimal"  formulas  have  been  obtained  which  mini- 
mize |R  (xP+l>|  among  all  symmetric  Chebyshev  quadratures  of  degree 
n 

p = 2[n/2]  - 1 . We  show  here  that  these  same  formulas  in  fact  minimize 

| R (x*)  | for  each  i > p + 1 . 
n — 
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ON  OPTIMAL  CHEBYSHEV-TYPE  QUADRATURES 


Walter  Gautschi  and  Giovanni  Monegato 


1.  Introduction.  We  consider  symmetric  Chebyshev-type  quadrature 
formulas  of  " next-to-highest"  algebraic  degree  of  exactness,  i.e.,  equally 
weighted  quadrature  formulas  of  the  form 

1 2 n 

(1.!)  / f(x)dx  = - Yj  f(tk>  + R(0  , 

-1  n k^l  k n 


subject  to  the  following  constraints, 


( all  t.  are  real,  t > t > . . . > t , 
i k 1 — 2 — — n ’ 


(1.2) 


/ fn+l-k  + *k  = °»  k = 1,  2,  . . . , n (symmetry), 


Rn(f)  = 0 for  every  polynomial  f of  degree  p = 2[n/2]  - 1 . 

If  1 < n < 7,  or  n = 9,  the  classical  Chebyshev  quadratures  satisfy  all 
these  conditions,  the  last  one  in  the  strengthened  form  with  p = 2[n/2]  + 1 . 
Other  Chebyshev-type  formulas,  satisfying  (1.  2),  exist  only  for  n = 8,  10, 
11,  and  13.  This  is  shown  in  [6]  and  [1],  where  formulas  are  derived  that 
are  "optimal"  in  the  sense  of  minimizing 


(1.  3) 


t)  = 


"v*p+1> 


p = 2[n/2]  - 1 , 
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subject  to  (1.  2).  Each  of  these  formulas  has  one  pair  of  symmetric  double 


1 


nodes,  or  one  double  node  at  the  origin,  which  (regretfully)  means  loss  of 
the  equal  coefficient  property.  We  will  show  here  that  these  same  formulas 
are  in  fact  optimal  in  the  much  wider  sense  of  minimizing 

(1.4)  Pi(t1,  t2,  . . . , tn)  = I Rn(>cl)  | for  each  i>p  + l , 

subject  to  (1.2).  (For  odd  integers  i,  the  statement  is  trivially  true, 
since  p = 0 by  virtue  of  symmetry. ) In  particular,  they  are  minimum  norm 
quadratures  in  the  Hardy  space  of  functions  f(z)  analytic  in  |z|  < r , 
r > 1,  and  square-integrable  on  | z | = r,  since  these  quadratures  minimize 

DO 

V r'2j[R  (xj)]2  (see,  e.  g. , [2]). 

J = P+1 

As  a preliminary  step  toward  these  results,  we  first  show  that  for 

each  n = 8,  10,  11,  13,  the  set  of  all  Chebyshev-type  quadrature  formulas 

(1.1),  satisfying  (1.  2),  forms  a one-parameter  family,  with  parameter  a in 

some  finite  closed  interval  X.  < a < . By  using  symbolic  Sturm  sequences, 

we  are  able  to  identify  X^  and  as  roots  of  certain  algebraic  equations. 

The  quadrature  formulas  corresponding  to  X < a < p all  have  n distinct 

n n 

(simple)  nodes  in  the  interval  (-1,1),  while  those  corresponding  to  a = X^ 
or  a = pn  have  a pair  of  symmetric  double  nodes,  or  a double  node  at  the 
origin.  The  main  result  then  follows  from  a (possibly  new)  monotonicity 
property  satisfied  by  the  power  sums  of  the  zeros  of  a polynomial.  In  par- 
ticular, it  transpires  that  the  optimal  formulas  [in  the  general  sense  of 
minimizing  (1.4)]  always  occur  at  a = X^  . They  are  necessarily  identical 
with  the  optimal  formulas  found  in  [6]  and  [1]. 
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2 . One-parameter  families  of  Chebyshev-type  quadratures.  We 


associate  with  (1. 1)  the  polynomial 


tt  n n-1 

|(x)  = I I (x-t  ) = x + a.x  + . . . + a 

k=l  k n 


From  Newton's  identities  it  follows  that  the  (real  ordered)  nodes  t^  will 
satisfy  (1.  2)  if  and  only  if 


r n n-2  2 

x + a^x  + . . . + a _x  + a n even)  , 
2 n-2  ’ 


(2.1)  |(x)=/ 


n n-2  3 . ... 

x + a2x  + . . . + a^  ^x  + ax  (n  odd)  , 


where  a is  some  real  parameter  and  the  coefficients  a_,  a^,  ....  a . 

2 4 p-1 

(p  = 2[n/2]  - 1)  are  given  uniquely  by 


s2  4 2a2  = 0 , 

S4  + a2S2  + 4a4  = ° ’ 


s , + a_s  + ...  + a s + (p-l)a  = 0 , 
p-1  2 p-3  p-3  2 p-1 


where 


= I / t2jdt  = 2?r,  * j = l,2,...,(P-l)/2  . 


2j  2 


-1 


2 j +1 


It  remains  to  determine  the  conditions  on  a under  which  the  polynomial  £ 


in  (2. 1)  has  only  real  zeros.  Letting 


£(*)=< 


’ * 2 

II  (X  ) (n  even)  , 


5{s  ^ 

(x  ) (n  odd)  , 
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we  are  seeking  the  conditions  on  a under  which  the  polynomial 


<2.  3) 


* U V - 1 

£ (x)  = x + a2x  + 


+ a2v  2x  + o,  v = [n/2]  , 


has  only  nonnegative  zeros.  While  the  problem  could  easily  be  discussed 
in  geometric  terms,  we  prefer  to  take  an  algebraic  approach  via  Sturm 
sequences.  Among  other  things,  this  allows  us  to  identify  the  limits  on 
the  parameter  a as  roots  of  certain  algebraic  equations  with  integer  co- 
efficients. 

We  recall  that  the  Sturm  sequence  of  a polynomial  pQ  of  degree  v is 
the  result  of  applying  Euclid's  algorithm  to  pQ  and  Pj  = pj^,  where  p^  is 
the  derivative  of  p^  . Thus, 


(2.4) 


P0  * Vi  ' P2  ’ 


P1  q2P2  ‘ P3  ’ 


PT-2  = qT-lPT-l  - PT 


^PT-l  = qrPT  » 


where  t < v,  and  each  p , 2 < k < t,  is  the  negative  remainder  of  the 
division  of  p , by  p . . Sturm's  theorem  (see,  e.  g. , [3,  p.  4481) 

K -2  K -1 

then  asserts  that  for  any  a < b,  if  PQ(a)  PQ(b)  * 0,  the  number  of  distinct 
zeros  of  pQ  in  [a,  b]  is  equal  to  v(a)  - v(b),  where  v(x),  for  fixed  x, 
is  the  number  of  sign  changes  in  the  numerical  sequence  PQ(x),  p^x),  . . ., 
Pt(x)  . Moreover,  if  xQ  is  a (real  or  complex)  zero  of  pQ,  then  its 
multiplicity  is  m if  and  only  if  xQ  is  a zero  of  the  terminal  polynomial 
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of  multiplicity  m-1  . In  particular,  all  zeros  of  pQ  in  [a,  b]  are 
simple  if  and  only  if  p ^ has  no  zeros  in  [a,b]  . 

* 

In  wishing  to  apply  this  to  the  polynomial  pQ  = £ in  (2.  3),  one  has 
to  be  prepared  to  operate  on  polynomials  whose  coefficients  depend  rationally 
on  the  parameter  a . This  is  best  done  with  the  help  of  symbolic  formula 
manipulation  systems.  We  indeed  generated  the  Sturm  sequences  symbolic- 
ally, using  the  MACSYMA  system  available  to  one  of  us  at  Stanford 
University,  and  repeating  the  computations  with  the  SAC-1  system  on  the 

UNIVAC  1110  computer  at  the  University  of  Wisconsin.  The  results  of  both 

| 

computations  were  identical,  and  are  summarized  below. 

Writing 

, . V-K  V-K-l  „ . 

p (x)  = c x + c , x + . . . + c . k = 0, 1,  . . . ,v  , 

K KK  K,K+1  Ky  > > > > 

the  coefficients  c , k < are  polynomials  or  rational  functions  of  a , 
with  rational  coefficients,  the  degrees  of  which  are  as  indicated  in  the 
schedules  below.  Zero  degree  always  represents  a non-vanishing  (rational) 
constant;  slashes  separate  numerator  degree  from  denominator  degree  in  case 
of  rational  functions.  In  all  rational  functions  on  the  diagonal,  the  denomi- 
nator polynomial  is  the  square  of  the  numerator  polynomial  in  the  rational 


n=8  ( v=4) 

0 0 0 0 1 
0 0 0 0 
0 0 1 
1 1 
3/2 


n=l0, 11  (v=5) 

0 0 0 0 0 1 

0 0 0 0 0 

0 0 0 1 

0 1 1 

2 2 
4/4 


n=l3  (v=6) 

0 0 0 0 0 0 1 

0 0 0 0 0 0 

0 0 0 0 1 

0 0 11 
1 1 1 

3/2  3/2 
7/6 


A 


function  preceding  it  on  the  diagonal.  The  numerator  polynomials  in  the  term- 
inal entries  are  irreducible  over  the  field  of  rationals,  except  in  the  case  of 

n = 13.  where  the  numerator  polynomial  of  c, , factors  into  the  square  of  a 

no 

linear  polynomial  and  a polynomial  of  degree  five. 

5{f 

According  to  Sturm's  theorem,  the  polynomial  £,  has  v distinct 
positive  zeros  (necessarily  all  simple)  if  and  only  if  v(0)  = v and  v(<*>)  = 0, 
that  is,  if  and  only  if 


(2.  5)  c > 0 and  (-1)V  * c >0  for  k = 0, 1,  2,  . . . , v . 

KK  K V 


By  examining  the  explicit  form  of  these  (polynomial)  inequalities  for  the 
four  values  of  n under  study,  one  finds  that  (2.  5)  is  equivalent  to 


(2.6) 


X.  < a-  < p , 
n n 


where  \ and  u are  either  zero,  or  equal  to  one  of  the  zeros 
n n ’ 

afn^  < < . . . of  c (a),  as  indicated  in  Table  2. 1.  Numerical 

12  3 vv 


8 

10 

11 

13 


0 

,(10) 

2 

4"’ 

(13) 


a 


(8) 

2 

0 

<n, 

(13) 


Table  2.1.  Parameter  intervals  for  symmetric 
Chebyshev-type  formulas  of  next-to-highest 
degree 
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N 

values  of  X , u , and  the  numerator  polynomial  c of  c in  exact 
n n vv  vv 

integer  form  (of  which  arf,n\  a^  are  zeros)  are  shown  below.  + 

By  verifying  that  v(l)  = v(oo)  = 0 for  a in  the  interval  (2.6),  one 

* 

finds  that  all  v zeros  of  £ are  in  fact  strictly  between  0 and  1 . 


n = 8 


CN  (°) 
44 


10766  1642834375a3  . 17  39295676  125a2  + 1639058085a  + 1698929 


= 0 


1.696315226301022019273508  X 10‘ 

n = 10 


<> 


10 


10 


N 


C55(a) 


103836101679541153125a4  - 699972937051972500a3 

- 118379410645986a2  + 59483935980a  + 10026277 

- 1.70502634454  1702458795446  X 10“4 
0 

n = 11 

4355201766189061767 168000000000000000a4 

-13748658484054726502645760000000000a3 

-32517908999837432625626726400000a2 

-17819369567349884423809920000a  - 2271280061895695934118607 

-4 


11 


11 


= -7.231729443377509440273227  X 10 


= -1.  84  3814  3654057  86592111158  X 10 


-4 


N 


t For  n = 13,  only  the  factor  (of  degree  5)  of  c is  shown  which  is 
relevant. 
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1 


n = 1 3 


# 


c.N,  (o)  = -3  1 194  312  382  3800644900  10  30897  982  132062886297600000000000000c5 

66  4 

+ 17  2 04  0 3029447  37  177  389  37  992  5 598908  3 36177295  36000000000000c 

-903866427  105482  168767  1987  188136769  191568998400000000c3 

7 

+ 1802928  58096  1992  1776960292  38  318998  36958  3104  000000c 
-15569598749 1956246845496029654386927806054400c 
+49094  135  59683  598477  099  364  3 50584692  8157  57 

X.  = 1.  02426094095810347493167  1 X 10'4 

1+13  = 1.  3259792845853337  15590553  X 10"4  . 

The  Sturm  sequence  terminates  prematurely  [i.e.  , t < v in  (2.4)] 
if  and  only  if  a is  a zero  of  one  of  the  diagonal  entries  c . We  ex- 

KK 

amined  the  Sturm  sequence  in  each  one  of  these  cases  and  found,  as  ex- 
pected, that  £ has  the  desired  number  of  nonnegative  zeros  precisely 
when  a is  a nonzero  endpoint  of  the  parameter  intervals  in  Table  2.  1. 

In  each  case,  £ has  one  double  zero  and  v - 2 simple  zeros  in  the 
open  interval  (0, 1)  . 

The  remaining  endpoint  a = 0 in  the  cases  n = 8 and  n = 10  can 

be  handled  by  forming  the  Sturm  sequences  at  x = « > 0 sufficiently  small, 

* 

and  at  x = 1 . One  finds  v - 1 simple  zeros  of  £ in  (0,1),  which 
together  with  the  zero  at  x = 0 again  yields  the  desired  number  of  zeros 
in  [0,1)  . 

We  may  summarize  our  findings  in  terms  of  the  polynomial  £ in 
(2.  1)  as  follows;  If_  n = 8,  10,  11,  or  13,  the  polynomial  £ has  n real 
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( symmetric)  zeros,  counting  multiplicities,  exactly  if  a is  in  the  interval 


\ < a < u . If  a is  an  interior  point  of  that  interval,  all  zeros  of  £ are 
n - — n — 1 — 

in  fact  simple  and  contained  in  the  open  interval  (-1,1).  I_f_  a is  one  of 

the  endpoints  of  f\  , p 1,  then  £ has  a pair  of  (symmetric)  double  zeros 

n n 

in  (-1,  1),  if  at  0,  and  a double  zero  at  the  origin,  if  a = 0,  all  other 

zeros  being  simple  and  located  in  (-1,  1)  . 

The  intervals  f V , p 1 found  here  are  in  agreement  with  the  (less 
1 n n 

accurate)  intervals  given  in  [4,  Table  I]  for  n = 8,  10,  and  11.  Our  results 
for  n =8  and  n = 10  contradict  a theorem  of  Pecka  [5],  which  indeed  is 
false  because  of  computational  errors  in  the  proof. 

3.  Optimal  Chebyshev-type  formulas  . In  each  one -parameter 
family  of  Chebyshev-type  quadratures,  obtained  in  Section  2,  we  now  wish 
to  single  out  one  that  minimizes  the  objective  function  in  (1.4).  We 
may  assume  i = 2j,  with  j > v,  where  v = [n/2]  . Letting  t (a)  de- 

’ K 

note  the  zeros  of  £ , and 

v 

sjV>  = Z £0c*)]; » j = o,i,2,...  , 

K =1 

their  power  sums,  and  assuming  \ <_  a <c  p.^,  the  objective  function  is 
easily  found  to  be 

(3.i)  = -iiaV1'  f- v • 

jjc 

We  thus  wish  to  minimize  p ^ ^ ( o)  on  the  interval  \n  < a < pn  . The 
lemmas  which  follow  will  lead  us  to  the  solution  of  this  problem. 


r 


Lemma  3. 1.  Suppose  the  algebraic  equation 


. . m m-1 

g(x)  :=  x + c,x  +...  + c ,x  + c =0 
1 m-1  m 


with  real  coefficients  c , has  m distinct  nonneaative  roots  x.  - x (c.. 

p - k k 1’ 


C2’  ' ' ' ’ Crr^’  k ' 1,  2,  . . . , m . Let  s denote  their  j -th  power  sum, 


m 


S)(Ci'C2 C*,\VVcl-V""cm>l'  ) = 0,1,2 


Then 


(3.2) 


9s  as 

— < 0 |f  j > p,  and  --7  = 0 if_  j < p . 


dc 


Proof.  Since,  by  Newton's  identities,  s.  depends  only  on 


cl>  C2>  • • • > Cj>  the  second  half  of  (3.  2)  is  self-evident.  It  suffices, 
therefore,  to  consider  j > p,  1 < p < m . 


Differentiating  the  identity  g(x,  (c.,  c_,  . . . , c ))  = 0 with  respect  to 

x i c m 


c , one  finds 
P 


3\  xk 


m-p 


m-p 


dc^  " ' g'(xJ 


k IT  (Xk  - X , 

l *k 


hence 


9s 


(3.  3) 


m . . 9x. 

J _ , V k . ri 

dc  k 9c  J /i 

P k=l  p 


m x 


m-p+j-1 


k=l  ,7'vv 

f *k 


The  sum  on  the  far  right  is  precisely  the  (m-l)-st  divided  difference 


[xl>  x2>  • • • > xm]f  of  the  function  f(x)  = x"1'^ -1  . Since 
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If  = f^m"^(x)/(m-l) ! = (m“^  *)  x^"^,  where  x is  in  the 
1’  2’  ’ mJ  m-1  ’ 


open  interval  spanned  by  x^,  x^,  . . . , x , and  since  all  x^  > 0,  it  follows 


that  x > 0,  and  the  sum  in  question  is  positive.  This  proves  the  lemma. 

3s, 


We  remark  that  the  proof  of  Lemma  3.  1 also  yields  >0  if  j < 0 , 


if  we  assume  all  x,  > 0,  but  we  will  not  need  this  fact  in  what  follows, 
k 


Lemma  3.  2.  If  n = 8,  10,  11,  or  13,  we  have 


(3.4) 


SAK)  < n 


jl'V  ~ 2(2J+1)  — j - V 


Proof.  Since  the  quadrature  formula  (1.  1)  has  polynomial  degree  of 


2v-2 


exactness  p = 2y-l,  it  is  certainly  exact  for  f(x)  = x , which  yields 


(3.  5) 


v-1  2(2v-l) 


For  any  j > v,  we  have 


# m # V-1  'e  1-V+l  * 

Sj  («)  = S [M«)]  [tK(«)l  <[t1(or)]J  SvA,  kn<a<Pn 

K=1 


# £ 

where  t^  (a)  is  the  largest  among  the  roots  t (a)  . The  lemma  will  hold 


for  any  value  of  j for  which 


r * .J-v+1  * . n 

vi* 


2(2j+l) 


By  virtue  of  (3.  5),  this  is  equivalent  to 


1 


# 2 v - 1 j - v+1 


.li. 


(3.6) 


An  elementary  computation  shows  that  the  right-hand  expression  in 

(3.6)  is  an  increasing  function  of  j for  j _>  v . By  consulting  the  numerical 
values  of  t (Xn>  in  [1],  and  noting  that  ^(X^)  = [t^X^]  , it  is  possible 

to  verify  the  validity  of  (3.  6)  for  j = v+2  when  n = 8,  for  j = v + 3 when 
n = lo,  for  ) = v + 4 when  n = 11,  and  for  j = v + 5 when  n = 13  . The 
inequality  therefore  holds  for  these  and  all  larger  values  of  j . For  the 
few  remaining  values  of  j,  one  can  verify  (3.4)  directly,  from  the  data 
available  in  [1].  This  completes  the  proof  of  Lemma  3.  2. 

The  problem  of  minimizing  P2.(a)  in  is  nOW  easily  solved- 

Theorem  3.  1.  If  n = 8,  10,  11,  or  13,  then  for  each  j > v , 

(3.7)  P^V  < P^")  — Xn<°- 

Proof.  Applying  Lemma  3. 1 to  the  polyncniial  | in  (2.  3),  the 

zeros  of  which  are  t*(  a),  we  find  that  each  s (a),  j > v,  is  strictly 

K J 

decreasing  on  \^<  a < pn,  hence 

(3.8)  s*(o)  < 8jV(Xn)  for  \n  < « < pn  . 

Combining  (3.  8)  with  Lemma  3.  2 gives 

for  each  j > v . Since  s*(a)  decreases  on  \^<  a < pn,  it  follows  that 

o*  la)  is  an  increasing  function  of  a,  hence  (3.7). 

K2j 

Corollary.  Let  {w. } be  any  sequence  of  weights  w^  > 0,  with 
> 0 for  some  j > v,  satisfying 
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^ j ~2  W < » , 


P " ( ° ) = Yj  wj[p2j(a)}2’  xn  i " 1 *V 

j = v 


P*(Xn>  < P*(a)  for  \ < ° < Pn  ‘ 


Proof.  Since  s*(  a)  < v[t*(a)]^>  where  ^(a)  is  the  largest  among 

the  roots  t%),  and  since  t'(a)<l  for  X <a<  p . the  sequence 
k 1 n ii 

. * . -1 
s.  (a)  tends  to  zero  at  least  geometrically,  and  consequently  p^lo)  ~ i 

as  j -*  oo  . The  infinite  series  in  (3.10),  therefore,  converges  by  virtue 

of  (3.  9),  and  (3.  11)  is  an  immediate  consequence  of  (3.  7). 
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